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摘 要 : 西藏 青 牺 主要 种 植 区 是 西藏 地 区 主要 粮食 作物 的 产 区 ,在 气候 变化 和 人 类 活动 加 剧 的 背景 
下 ,该 区 的 农业 深 受 干旱 影响 。 根 据 《 气 象 干旱 等 级 》(GB/T 20481-2017) 国 家 标准 ,计算 西藏 青 牺 


主要 种 植 区 25 个 气象 站 点 全 生育 期 的 逐日 气象 干 


旱 综合 指数 (MCI) ,统计 1981 一 2020 年 各 站 点 的 


气象 干旱 过 程 , 并 进一步 分 析 了 西藏 青 释 主要 种 植 区 各 生育 期 气象 干旱 的 时 空 变 化 特征 及 致 灾 危 
险 性 等 级 分 布 。 结 果 表 明 :(1) 近 40a 来 西藏 青 牺 主 要 种 植 区 呈现 气温 显著 上 升 (P<0.01) ,降水 不 
显著 增多 的 特点 ;季节 性 的 气温 和 降水 变化 趋势 差异 明显 。(2) 于 旱 覆 羡 范 围 和 持续 日 数 总 体 呈 现 


不 同 程度 的 下 降 趋 势 ,于 旱 强 度 ( 除 播种 -分 药 期 外 ) 均 有 所 增强 ;分 药 - 抽 入 期 干 时 最/ 


重 ,1981 一 


1990 年 干旱 最 重 ,2001 一 2010 年 干旱 最 轻 ,(3) 各 生育 期 干旱 发 生 频 率 的 高 值 区 主要 集中 在 日 路 
则 市 东部 、 拉 萨 市 南部 、 昌 都 市 西南 部 和 山南 市 中 部 地 区 , 低 值 区 位 于 林 艺 市 中 部 。(4) AER 
数 法 (EOF) 分 析 的 各 生育 期 主要 模 态 有 中 部 与 中 西部 相反 、 全 区 一 致 .西南 与 东北 相反 。(S$) 各 生 
育 期 干旱 致 灾 危 险 性 等 级 总 体 呈 现 西 高 东 低 的 分 布 特征 ,其 中 ,分 药 - 抽 穗 期 致 灾 危 险 性 总 体 较 


高 , 抽 德 -成 熟 期 次 之 ,播种 -分 药 期 最 低 。 


X 键 词 : MCI 指 数 ; 致 灾 危险 性 ; 干旱 强度 ; BR; 西藏 


文章 编号 : 1000 - 6060(2023)07 - 1098 -= 13( 1098 ~ 1110) 


随 着 全 球 变 暖 ,干旱 发 生 频 率 增多 .强度 增加 ， 
影响 范围 随 之 扩大 "1, 对 农 牧 业 、 生 态 环 境 水 资源 
和 社会 经 济 的 发 展 造成 了 不 可 估量 的 损失 。 因 此 ， 
干旱 指标 的 选取 .干旱 过 程 的 判断 .干旱 的 演变 规 
律 以 及 干旱 灾害 的 评估 等 问题 越 来 越 受 到 国内 外 
专家 和 学 者 的 关注 和 重视 ”。 

根据 干旱 成 因 和 影响 ,将 干旱 分 为 气象 干旱 、 
TK SCF RAE BR AREAS SEY 而 气 
象 干旱 是 其 他 类 型 干旱 的 起 因 和 监测 评估 的 基 
础 。 目 前 ,国内 外 向 用 的 气象 干旱 指数 有 标准 化 降 
水 指数 (SPI)” \、 帕 默 尔 干旱 强度 指数 (PDSI)"\ 标 准 
化 降水 蒸 散 指数 (SPEI)5.Z 指数 “综合 气象 干旱 
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指数 (CU7” 等 。SPI 指 数 适用 于 所 有 气候 条 件 ,但 未 
考虑 气温 对 干旱 的 影响 ;PDSI 指 数 和 SPEI 指 数 均 考 
虑 了 降水 和 和 气温 ,缺点 是 都 需要 连续 完整 的 数据 ， 
时 间 尺 度 会 导致 干 旱 监测 滞后 ;Z 指 数 来 源 于 对 PD- 
SI 指数 的 改进 ,经 常用 于 不 同时 段 干 旱 的 比较 ,也 
用 来 计算 干旱 结束 的 时 间 , 但 在 中 国 使 用 较 少 ;CI 
指数 主要 计算 干旱 频率 和 持续 时 间 , 涉 及 数据 易 获 
取 , 但 是 在 近 几 年 重大 干旱 事件 的 监测 业务 中 暴露 
出 干旱 发 展 过 程 中 等 级 异常 跳 变 .严重 干旱 事件 反 
应 偏 轻 等 问题 。 基 于 此 ,本 文选 用 国家 气候 中 心 
推出 的 气象 干旱 综合 指数 (MCI) , 既 弥 补 了 CI 指数 
的 缺点 ,同时 考虑 到 近期 降水 和 长 期 降水 的 影响 、 
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冬 半 年 和 夏 半年 的 不 同 北方 和 南方 的 不 同 、 农 作 
物 的 生长 情况 以 及 当地 水 资源 的 情况 。 

西藏 气候 独特 且 复 杂 多 样 。 海 拔高 ,气温 低 ， 
昼夜 温差 大 ;降水 少 ,地 区 差异 大 ;太阳 光照 强 ,日 
照 时 间 长 ,是 导致 其 易 发 生 干 旱 的 主要 原因 。 据 统 
计 , 西 藏 1983 一 2013 年 共 发 生 旱灾 356 次 ,平均 每 
年 11 次 ,干旱 灾害 的 高 发 期 出 现在 春 夏季 ,也 是 青 
牺 的 主要 发 育 期 ”。 青 梨 是 西藏 主要 的 粮食 作物 ， 
主要 分 布 在 西藏 高 原 雅 鲁 藏 布 江 中 游 河 谷 的 拉萨 、 
日 喀 则 .山南 市 以 及 藏 东 三 江 流域 的 昌都 市 和 藏 东 
南 的 林芝 市 ,该 地 区 年 平均 气温 显著 升 高 ,而 降水 
总 量 微弱 增加 ,干旱 随 年 份 的 递增 有 加 重 的 趋势 "" 。 
青 穆 受到 的 干旱 影响 主要 来 自 降 水 量 和 蒸 散 发 量 ， 
相关 人 研究 表明 ,自然 降水 是 左右 其 旱情 的 决定 性 因 
素 ", 青 牺 在 分 更 .拔节 和 灌浆 期 耗 水 强度 最 大 , 降 
水 偏 少 会 影响 其 正常 生长 发 育 ,降低 青 梨 的 产量 与 
品质 “”。 对 比 其 他 农作物 的 情况 , 邓 国 卫 等 ” 认 
为 从 移 栽 到 成 熟 , 水 称 干 旱 频 率 随 着 生育 进程 未 渐 
降低 ,各 生育 期 干旱 对 水 称 产 量 的 影响 逐渐 减弱 ; 
黄岩 等 史 认 为 拔节 期 同等 程度 干旱 造成 的 玉米 减 
产 幅度 均 大 于 苗 期 和 灌浆 期 干旱 ; 曹 永 强 等 认为 
河北 省 夏 玉 米 干 旱 主要 发 生 在 拔节 -抽雄 和 抽雄 - 
成 熟 期 ; 吴 蓝 芒 ” 认 为 冬 前 ,冬小麦 干旱 发 生 频 率 
整体 较 低 ,而 冬 后 ,干旱 频率 随 生 育 期 的 推移 有 所 
WME, 可 见 不 同 作物 在 不 同 生育 期 发 生 的 干旱 概率 
和 减产 程度 是 不 尽 相 同 的 。 因 此 ,文中 考虑 了 青 牺 


高 程 m 
OER 2734.51 0 


不 同 生 育 期 的 干旱 演变 规律 。 

目前 ,针对 西藏 地 区 的 干旱 研究 大 多 停留 在 单 
项 气象 干旱 指数 的 干旱 演变 特征 分 析 ” ,时 间 尺 
度 过 于 单一 ,无 法 全 面 反 映 该 地 区 的 干旱 变化 状 
况 。 与 此 同时 ,基于 MCI 指 数 作为 西藏 地 区 干旱 致 
灾 因 子 的 评估 鲜 见 报道 。 本 文 利用 MCI 指 数 ,从 青 
牺 不 同 生育 期 角度 出 发 ,探讨 西藏 青 穆 主要 种 植 区 
1981 一 2020 年 气象 干旱 时 空 演变 规律 ,并 进一步 定 
义 干 旱 致 灾 和 危险 性 指数 ,掌握 该 地 区 干旱 致 灾 危 险 
性 的 分 布 特点 ,以 期 识别 和 评估 干旱 风险 ,为 政府 
防 灾 减 灾 和 农业 可 持续 发 展 提供 科学 支撑 。 


1 数据 与 方法 


1.1 数据 来 源 

根据 史 继 清 等 中 和 杜 军 等 3 对 西藏 青 牺 的 研 
究 ,选择 西藏 全 区 (阿里 地 区 和 那曲 市 除外 ) 青 牺 播 
种 面积 大 于 1000 hm 的 县 (区 ) 作 为 西藏 青 穆 主 要 种 
植 区 (图 1)。 西 藏青 牺 主要 种 植 区 25 个 站 点 的 气象 
资料 (1981 一 2020 年 ) 来 源 于 西藏 自治 区 气象 局 , 青 
PRE BEE Fy DY RP SP BEY oP BEF tT 
成 熟 期 ) 资 料 来 自 西 藏 自治 区 气候 中 心 生态 与 农业 
气象 研究 室 ( 表 1)。 
1.2 研究 方法 
1.2.1 气象 干旱 综合 指数 (MCI) 干旱 是 由 于 降水 
长 期 亏 缺 和 近期 亏 缺 综合 效应 累加 的 结果 ,MCI 指 


180 km 


注 :该 图 基于 西藏 自然 资源 厅 标 准 地 图 服务 网 站 下 载 的 审 图 号 为 藏 S2022)004 号 的 标准 地 图 制作 , 底 图 边界 无 修改 。 下 同 。 
图 1 西藏 译 称 主要 种 植 区 的 气象 站 点 分 布 


Fig. 1 Distribution of meteorological stations in highland barley growing areas of Tibet 
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表 1 西藏 青 穆 主要 种 植 区 生育 期 资料 
Tab. 1 Growth period data of highland barley growing areas in Tibet 
站 后 FeO BE TP eS 抽穗 -成 熟 期 播种 -成 熟 其 
Hig 5 月 中 旬 一 6 月 上 旬 ”6 月 上 旬 一 7 月 中 旬 ”7 月 中 旬 一 9 月 上 旬 ”5 月 中 旬 一 9 月 上 旬 


EH EMAR 

MBE ETR JER 

日 喀 则 YTE BARR 

林芝 、 米 林 、 波 密 、 察 隅 

加 查 HE .隆子 . 错 那 . 浪 卡 子 

TH BAZAT ER AE RS A 


4 月 中 旬 一 6 月 下 旬 
4 月 中 旬 一 5 月 中 旬 
4 月 下 旬 一 5 月 中 旬 
3 月 中 旬 一 4 月 下 旬 
4 月 上 旬 一 5 月 中 多 
4 月 上 旬 一 5 月 中 多 


数 考虑 了 60 d 内 的 有 效 降水 (权重 累积 降水 ) .30 d 
内 蒸 散 (相对 湿润 度 ) 以 及 季度 尺度 (90 d) 降 水 和 近 
半年 尺度 (150 qd) 降水 的 综合 影响 。MCI 指 数 计算 
AAMT: 

MCI = Ka X (a X SPIW s + b X MI,, +c X SPI, + w 

d x SPI <0) 

式 中 : Ka NED DAT ARAL, Xx HA 1—3 H A 11—12 
月 (0)、4 月 (0.6)、5 月 和 8 月 (1.0)、6 一 7 月 (1.2)、9 
月 (0.9)、10 月 (0.4); SPIW。 为 近 60 d 标 准 化 权重 
降水 系数 ; ML 为 近 30 d 相 对 湿润 度 指 数 ; SPI, 、 
SPI s 分 别 为 近 90 d 和 近 150 d 标 准 化 降水 指数 ; a 
b c d IHA SPIW o ~ Mlo ~ SPloo ~ SPL so 的 权重 系 
数 , 这 里 取 值 分 别 为 0.3、0.5、0.3、0.2。 根 据 MCI 指 
数 划 分 的 气象 干旱 等 级 见 表 2。 
1.2.2 干旱 覆盖 范围 .发 生 频 率 和 干旱 强度 “o 
究 时 段 内 至 少 出 现 一 次 干旱 过 程 , 旦 累计 评价 时 段 
内 干旱 持续 时 间 超 过 所 评价 时 段 的 1/4 时 , 则 认为 
该 时 段 发 生 干 旱 事件 ”。 当 MCI 指 数 至 少 连续 10 d 
为 轻 旱 及 以 上 等 级 , 则 确定 为 一 次 干旱 过 程 。 干 旱 
过 程 时 段 内 第 一 次 出 现 轻 旱 的 日 期 ,为 干旱 过 程 开 
始 日 ; 当 MCI 指 数 连续 10 d 判 别 为 无 旱 等 级 时 , 则 
干旱 过 程 结束 , 将 最 后 一 次 MCI 指 数 达 到 无 早 等 级 
时 的 日 期 定义 为 干旱 过 程 结束 日 期 。 从 干旱 过 程 
开始 到 结束 期 间 所 持续 的 时 间 为 干旱 持续 时 间 。 
将 一 次 干旱 过 程 内 累计 的 轻 旱 及 以 上 等 级 MCI 指 
数 之 和 定义 为 这 次 干旱 过 程 的 强度 。 

基于 站 点 计算 干旱 履 盖 范围 ($) ,计算 公式 


表 2 气象 干旱 综合 


6 月 下 旬 一 7 月 下 旬 
5 月 中 旬 一 6 月 中 多 
5 月 中 旬 一 6 月 下 旬 
4 月 下 名 一 6 月 中 名 
5 月 中 旬 一 6 月 下 旬 
5 月 中 旬 一 7 月 上 名 


7 月 下 旬 一 9 月 中 旬 
6 月 中 旬 一 8 月 上 外 
6 月 下 旬 一 8 月 中 旬 
6 月 中 旬 一 7 月 下 旬 
6 月 下 旬 一 8 月 上 旬 
7 月 上 旬 一 8 月 中 旬 


4 月 中 旬 一 9 月 中 旬 
4 月 中 旬 一 8 月 上 名 
4 月 下 旬 一 8 月 中 名 
3 月 中 旬 一 7 月 下 旬 
4 月 上 旬 一 8 月 上 名 
4 月 上 名 一 8 月 中 名 


如 下 : 
S= y x 100% (2) 
AP: m 为 每 年 发 生 干 旱 事件 的 站 点 个 数 ; WM 为 气 
象 站 点 总 个 数 ,MM 取 值 25。 
干旱 发 生 频 率 (P) 计 算 公 式 如 下 : 
P=- x 100% (3) 
SUH: n 为 发 生 干 旱 的 年 数 ; N 为 研究 年 份 , N 取 值 
40。 
干旱 强度 (D) 计 算 公 式 如 下 : 
DO=Y1 (4) 


式 中 :为 一 次 干旱 过 程 内 的 干旱 天 数 ; /为 第 i 天 
MCI FEX, SiS ko 

1.23 干旱 致 灾 危 险 性 干旱 致 灾 和 危险 性 是 由 旱灾 
的 强度 和 频率 共同 决定 的 ,本 文 基 于 干旱 发 生 频 率 
和 干旱 强度 定义 干旱 致 灾 危 险 性 指数 (R) ,其 计算 
AAJ: 


R=PxD (5) 
式 中 : P 为 各 站 点 的 干旱 发 生 频 率 ; D 为 各 站 点 的 
干旱 强度 。 
1.2.4 其 他 方法 采用 气候 倾向 率 法 人 研究 气候 变化 
和 MCI 指 数 的 时 间 变 化 趋势 ;选用 Kriging 插 值 对 干 
旱 强 度 和 和 干旱 致 灾 危 险 性 进行 空间 插值 ;对 干旱 强 
度 进行 经 验 正 交 函数 法 (EOF) 分 析 ”。 
EOF 的 分 解 原理 是 将 有 限 空间 内 呈 不 规则 分 
布 的 观测 站 点 数据 ,分 解 为 少数 互相 独立 的 典型 模 


EA (MCT) 等 级 划分 


Tab.2 Classification of drought with meteorological drought composite index (MCI) 


干旱 等 级 无 旱 轻 旱 
MCI 取 值 MCI>-0.5 -1.0<MCI<-0.5 


-1.5<MCI<-1.0 


中 早 重 早 om 


-2.0<MCI<-1.5 MCI<-2.0 
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态 ,该 方法 既 可 以 取 空间 不 同 站 点 进行 分 解 ,也 可 
以 对 同一 站 点 的 不 同时 间 ,不同 高 度 的 多 种 要 素 进 
行 综 和 分 析 , 现 广泛 应 用 于 气象 学 领域 "; 同 时 ， 
EOF 分 解 的 展开 收敛 速度 均 很 快 ,变量 场 的 信息 
能 够 很 容易 地 集中 在 前 几 个 模 态 上 ;累计 方差 贡献 
率 达 到 一 定 程度 的 前 几 个 模 态 具有 一 定 物 理 意义 ， 
可 有 效 反 映 干旱 强度 时 空 场 的 变化 特征 ”1, 基 本 思 
路 是 将 包含 g 个 空间 点 个 年 份 的 观测 值 构成 的 要 素 
AOX a ) 分 解 为 相互 正 交 的 空间 函数 矩阵 (了 ) 
FINTE RUER CT, ) 的 乘积 , 即 : 


g 


A eh = YD Vex, x Ta) p= l, 2, & (6) 


p=1 


其 中 ， 
Vi Vig OU 
re 3 T (7) 
Pala ° Vap 
taito e bp 
Tai = e e R (8) 
tatn e thy 


2 结果 与 分 析 


21 西藏 青 牺 主 要 种 植 区 气温 、 降 水 变化 趋势 

由 图 2 可 知 ,1981 一 2020 年 西藏 青 穆 主要 种 植 
区 多 年 平均 降水 量 为 490.92 mm, 上 且 以 7.71 mm- (10a) ' 
的 速率 呈 不 显著 增加 趋势 ,其 中 ,1983 年 和 1998 年 分 
别 是 近 40 a 该 研究 区 年 平均 降水 量 最 少 ( 350.11 mm) 
和 最 多 (611.20 mm ) 的 年 份 。 该 研究 区 多 年 平均 气 
温 为 6.52 € ,并 以 0.34 C- (10a) HEK E m E E 


CO 年 平均 气温 
750 NETE DEET 线性 (年 平均 气温 ) 
Jia-0.0341xz+5.8246 
R?=0.5783 


年 平均 降水 量 /mm 


升 趋势 (P<0.01) ,其 中 1997 年 和 2009 年 分 别 是 过 去 
40 a 该 研究 区 年 平均 气温 最 低 (5.46 °C ) 和 最 高 
(7.58 C) AES 

由 图 3 可 知 , 近 40a 西 藏青 牺 主 要 种 植 区 春季 
多 年 平均 降水 量 为 84.99 mm, 占 年 平均 降水 量 的 
17.31% ,并 以 4.83 mm:(10a)' 的 速率 呈 显 著 增加 趋 
势 (P<0.05) ;夏季 多 年 平均 降水 量 为 296.86 mm , 占 
年 平均 降水 量 的 60.47%, 并 以 6.71 mm + (10a)! HJ ZE 
率 呈 不 显著 增加 趋势 ; 秋 、 冬 季 多 年 平均 降水 量 分 
别 为 94.69 mm 和 14.39 mm ,分 别 以 -3.22 mm: (10a)! 
和 -0.42 mm (10a) "的 速率 呈 不 显著 的 减少 趋势 。 
春 、 夏 、 秋 、 冬 季 多 年 平均 气温 分 别 为 6.58 C. 
14.04% .7.05% 和 -1.63 % ,并 分 别 以 0.27 (10a)'、 
0.26 C- (10a) .0.38%.(10a 0.46 C- (10a) HEX 
呈现 显著 的 上 升 趋势 (P<0.01)。 
2.2 干旱 覆盖 范围 的 年 际 变化 特征 

1981 一 2020 年 西藏 青 穆 主 要 种 植 区 各 生育 期 
干旱 上 覆盖 范围 的 变化 趋势 如 图 4 所 示 。 

播种 -分 葛 期 和 抽穗 -成 熟 期 ,干旱 覆盖 范围 均 
呈 不 明显 的 下 降 趋 势 ,多 年 均值 分 别 为 42.7% 和 
31.4%。 播 种 -分 堤 期 ,各 年 份 均 有 一 定 范 围 的 干旱 
发 生 ; 从 各 年 代 来 看 ,1991 一 2000 年 平均 干旱 覆盖 
范围 最 广 ,为 51.6% ,200] 一 2010 年 最 小 ,为 32.4% 。 
抽穗 -成 熟 期 ,1983 年 干旱 履 善 范围 最 广 ,为 84.0%， 
2000, 2004 4E V4 Æ 2016—2017 年 干旱 覆盖 范围 为 
0, 表明 这 些 年 份 研究 区 无 干旱 发 生 ; 各 年 代 中 ， 
1981 一 1990 年 干旱 发 生 范 围 最 大 ,为 43.3% ,2001 一 
2010 年 最 小 ,为 24.0%。 


一 一 一 年 平均 降水 量 

aaa ol: 线性 (年 平均 降水 量 ) 10 
Vee =0.7712x+475.11 

R°=0.0192 


EEC 


0 0 
1981 1984 1987 1990 1993 1996 1999 2002 2005 2008 2011 2014 2017 2020 
年 份 


图 2 1981 一 2020 年 西藏 青 牺 主 要 种 植 区 年 平均 降水 量 与 年 平均 气温 的 变化 
Fig. 2 Changes of annual mean precipitation and annual mean temperature 


in highland barley growing areas of Tibet from 1981 to 2020 
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图 3 1981 一 2020 年 西藏 青 称 主 要 种 植 区 季节 平均 降水 量 与 平均 气温 的 变化 
Fig. 3 Changes of seasonal mean precipitation and mean temperature 


in highland barley growing areas of Tibet from 1981 to 2020 
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图 4 1981 一 2020 年 西藏 青 牺 主 要 种 植 区 各 生育 期 干旱 履 盖 范 围 的 变化 
Fig. 4 Change of drought coverage in different growth stages in highland barley growing areas of Tibet from 1981 to 2020 
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Oy Be — FH A] A ee W, A aoe FY 
现 显著 (P<0.05) 的 下 降 趋 势 ,多 年 均值 分 别 为 
51.3% 和 36.8%。 分 也 -抽穗 期 共有 6 af my0 Al 
过 80.0% ,其 中 ,1982 年 和 1986 年 覆盖 范围 均 为 
88.0% ,发 生 范围 最 广 ,2000 年 无 干旱 发 生 ; 各 年 代 
中 , 1981—1990 年 干旱 最 普遍 ,平均 上 覆盖 范围 为 
72.0% , 2001—2010 年 平均 履 盖 范围 最 少 , 为 
32.4% 。 全 生育 期 ,1992 年 干旱 覆盖 范围 最 广 ,为 
92.0% , 2004 4F Fl 2008 年 无 干旱 发 生 ,其 余年 份 干 
早 禾 盖 疙 围 不 一 ;各 年 代 中 ,1981 一 1990 年 干旱 履 
盖 范 围 最 广 , 为 54.2% , 2001 一 2010 年 最 小 , 仅 为 
19.6%. 

2.3 干旱 持续 日 数 的 年 际 变化 特征 

由 图 5 可 知 ,播种 -分 碍 期 .分 蓝 - 抽 穗 期 抽穗- 
成 熟 期 和 全 生育 期 的 干旱 持续 日 数 均 呈现 不 显著 
的 下 降 趋 势 ,多 年 均值 分 别 为 24d.26 d、25 d 和 60 d- 

播种 -分 葵 期 2000 年 干旱 持续 日 数 最 长 ,为 
35 d,2004 年 最 短 ,为 14 d; 各 年 代 干 旱 持 续 日 数 介 
于 23~25 d 之 间 。 分 幕 - 抽 穗 期 ,干旱 持续 日 数 最 长 
的 是 1988 年 的 39 d ,最短 的 是 2000 年 的 0 d, AFF 
发 生 ; 各 年 代 中 ,1981 一 1990 年 干旱 平均 持续 日 数 


一 一 干旱 持续 日 数 


(a) RPS EE 
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最 长 , 约 30 d,2001 一 2010 年 最 短 , 约 20 d, 其 余年 代 
均 为 28 d。 抽 穗 -成 熟 期 1983 年 干旱 持续 日 数 最 
长 ,为 49 d,2000.2004 年 和 2016 一 2017 年 均 无 干旱 
发 生 ; 各 年 代 中 ,1981 一 1990 年 干旱 平均 持续 日 数 
最 长 , 约 31 d,2001 一 2010 年 最 短 , 约 21 d。 全 生育 
期 ,干旱 持续 日 数 最 长 的 是 1984 年 的 136 d, 最 短 的 
是 2000 年 和 2008 年 的 0 d; 各 年 代 中 ,1981 一 1990 
年 干旱 平均 持续 日 数 最 长 , 约 73 d, 2001—2010 年 
最 短 , 约 44 d, 其 余年 代 在 55 d 左 右 。 
2.4 干旱 频率 的 空间 分 布 特征 

由 图 6 可 以 看 出 ,整体 上 ,分 尧 -抽穗 期 干旱 频 
率 最 高 ,播种 -分 刻 期 次 之 ,抽穗 -成 熟 期 最 低 。 播 
种 -分 星期 ,林芝 市 、 错 那 站 的 干旱 发 生 率 普 忆 偶 
低 , 均 在 30.0% 以 下 ;研究 区 东部 .西部 干旱 发 生 率 
多 在 40.0%~60.0% 之 间 , 占 总 站 点 数 的 68.0% 。 分 
Be — HH AY] , JP $F 50.0%~70.0% 之 间 的 站 点 数 最 多 ， 
占 比 达 64.0% , 主要 分 布 在 研究 区 中 西部 和 昌都 市 
中 南部 ;干旱 发 生 率 在 30.0% 以 下 和 70.0% 以 上 的 站 
点 均 为 1 个 ,分 别 出 现在 波 密 站 和 江 孜 站 。 抽 穗 - 成 
熟 期 ,干旱 发 生 率 均 在 50.0% 以 下 ( 除 隆 子 站 外 ) , 低 
值 区 主要 位 于 昌都 市 大 部 、 林 芝 市 中 部 地 区 。 全 生 


EE 线性 (干旱 持续 日 数 ) 
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图 5 1981 一 2020 年 西藏 青 牺 主 要 种 植 区 各 生育 期 的 干旱 持续 日 数 的 变化 趋势 
Fig. 5 Change trends of drought duration days in different growth stages of highland barley 


growing areas in Tibet from 1981 to 2020 
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图 6 1981—2020 E PIRE PR SE Be PAL EX 5 A IP 8 F928 A 
Fig. 6 Spatial distributions of drought frequency in different growth stages of highland barley 


growing areas in Tibet from 1981 to 2020 


育 期 ,干旱 发 生 率 高 值 区 集中 在 日 喀 则 市 中 东部 、 
昌都 市 西部 以 及 山南 市 北部 。 
2.5 干旱 强度 的 时 空 特征 

FH EOF 对 西藏 青 穆 主要 种 植 区 各 生育 期 干 
旱 强 度 的 标准 值 进行 分 解 ( 表 3、 图 7) ,分 析 结 果 表 
HH ,播种 -分 更 期 .分 药 - 抽 和 初期 .抽穗 -成 熟 期 和 全 
生育 期 前 5 个 特征 向 量 特征 值 的 累计 贡献 率 分 别 达 
到 71.2% 72.1% .83.5% 和 79.3% , 其 中 播种 -分 坦 期 
和 全 生育 期 仅 有 第 一 个 特征 根 的 误差 范围 不 重叠 ， 
分 艺 - 抽 穗 期 和 抽穗 -成 熟 期 则 有 前 2 个 特征 根 的 误 
BBA. It North 显著 性 检验 ,播种 -分 龊 
期 (第 一 载 集 量 ) od BEF E , ta Bt) 
抽穗 -成 熟 期 (第 一 二 载荷 量 ) 和 全 生育 期 (第 一 载 
傈 量 ) 的 累计 贡献 率 分 别 接 近 23.1%、51.8%、68.7% 
和 35.8% ,可 以 很 好 地 解释 近 40 a 来 西藏 青 穆 主 要 
种 植 区 各 生育 期 干旱 强度 的 主要 时 空 分 布 特征 。 

播种 -分 药 期 (图 7a) ,第 一 模 态 (EOF1) 的 干 早 
强度 负 值 区 集中 在 雅鲁藏布江 一 线 大 部 林芝 站 、 
米 林 站 ,其 余 为 正 值 区 。 从 EOF1 时间 序 列 ( 图 8a) 
可 以 看 出 ,1989 一 1998 年 和 2007 一 2020 年 干旱 强度 
时 间 系 数 在 正 负 间 交 其 变化 ,说 明 这 2 个 时 间 段 研 


究 区 干旱 强度 的 强 弱 变化 明显 ;1998 一 2006 年 均 为 
正 值 ,表明 该 时 间 段 研究 区 干旱 强度 减弱 ,结合 图 2 
发 现 ,1998 一 2005 年 降水 量 偏 多 ,因此 对 干旱 强度 
产生 一 定 的 影响 。 

分 艺 - 抽 穗 期 .抽穗 -成 熟 期 和 全 生育 期 (网 7b、 
d\f) 的 EOF1 在 整个 研究 区 大 体 上 为 正 值 ( 除 抽穗 - 
成 熟 期 的 肢 拉 木 站 和 全 生育 期 的 丁 青 站 以 外 ), 表 
明 近 40 a 来 研究 区 在 这 3 个 生育 期 的 干旱 强度 变化 
趋势 空间 上 具有 高 度 的 一 致 性 , 即 干旱 强度 在 研究 
区 普遍 增强 或 减弱 。 高 值 中 心 主 要 位 于 雅 鲁 藏 布 
江 一 线 的 山南 市 北部 .拉萨 市 南部 和 日 喀 则 市 东 
部 ,是 干旱 强度 变化 反应 最 敏感 的 3 个 区 域 。 从 图 
8b~d 可 以 得 出 ,分 蓝 - 抽 穗 期 .抽穗 -成 熟 期 和 全 生 
育 期 的 时 间 系 数 分 别 在 1999 一 2004 年 .1995 一 2004 
年 和 1996 一 2004 年 为 正 值 , 表 明 这 3 个 时 间 段 研究 
区 干旱 强度 减弱 ,与 相同 时 段 的 降水 量 偏 多 有 一 定 
的 关系 ;其 余 时 间 段 时 间 系 数 正 负 交替 变化 。 

分 药 - 抽 穗 期 的 第 二 模 态 (EOF2)( 图 7c) 的 干 
旱 强 度 整 体 表 现 为 研究 区 中 部 正 值 中 心 与 东部 、 西 
部 边缘 负 值 中 心 呈 相反 的 空间 分 布 特征 ;山南 市 、 
日 喀 则 市 东部 拉萨 市 南部 、 林 芝 市 局 部 边缘 地 区 
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表 3 EOF 前 $ 个 载 和 荷 量 方差 贡献 率 及 累计 方差 贡献 率 


Tab. 3 Variance contribution rate and cumulative variance contribution rate of the first five loads of EOF 


生育 期 EOF HAR RHN 
献 率 /% 贡献 率 /% 
播种 -分 蓝 期 第 一 载荷 量 23.1 23.1 
oh tay E 18.8 41.9 
of = ay at 13.6 55.8 
第 四 载荷 量 8.6 64.4 
第 五 载荷 量 6.8 71.2 
抽穗 -成 熟 期 第 一 载荷 量 61.5 61.5 
PRIT E 7.2 68.7 
SP ERIT E 5.4 74.1 
第 四 载荷 量 4.8 78.9 
第 五 载荷 量 4.6 83.5 


注 :EOF 为 经 验 正 交 函数 法 。 


以 及 丁 青 站 是 正 值 区 ,其 余 为 负 值 区 。 结 合 图 8b 发 
现 ,2000 一 2006 年 和 2017 一 2020 年 主要 表现 为 研究 
区 中 部 干旱 强度 减弱 ,而 东部 、 西 部 边缘 增 强 ; 
2007 一 2012 年 呈现 相反 的 变化 趋势 。 

抽穗 -成 熟 期 的 第 二 模 态 (EOF2) (Al 7e) ky F 
旱 强 度 表现 为 “西南 负 、 东 北 正 ” 的 反 癌 分布 型 , 即 
西南 东北 差异 型 ; 负 值 区 主要 集中 在 雅鲁藏布江 一 
线 南 段 . 沿 喜马拉雅 山 一 带 以 及 左 贡 站 、 芒 康 站 、 察 
阳 站 ,其 他 属于 正 值 区 。 由 图 8c 发 现 ,1981 一 2020 
年 的 时 间 系 数 在 正 负 间 振 荡 ,表明 研究 区 干旱 强度 
年 际 间 的 变化 趋势 显著 。 
2.6 干旱 致 灾 危 险 性 等 级 的 空间 分 布 

All FA ArcGIS 10.4 软件 ,分 别 将 播种 -分 刻 期 分 
Be — Fy .抽穗 -成 熟 期 和 全 生育 期 的 干旱 致 灾 危 
险 性 划分 为 低 、 较 低 .中 等 、 较 高 .高 等 5 个 等 级 (图 
9)。 王 旱 致 灾 危 险 性 在 播种 -分 碍 期 ,大 体 上 沿 洛 
隆 一 八 宿 一 左 贡 一 线 ,西部 地 区 呈现 由 东 向 西 逐 渐 
加 重 的 趋势 ;分 万 -抽穗 期 旱灾 危险 性 均 偏 高 ; 抽 
穗 - 成 熟 期 ,整体 上 呈现 研究 区 西高 东 低 的 分 布 特 
征 ; 全 生育 期 旱灾 危险 性 处 于 中 等 或 偏 低 。 播 种 - 
分 药 期 ,干旱 致 灾 危 险 性 中 高 等 级 集中 分 布 在 日 喀 
则 市 中 部 ,其 他 地 区 危险 性 等 级 偏 低 , 其 中 林芝 市 
大 部 和 山南 市 南部 最 低 。 分 药 - 抽 穗 期 , 雅 鲁 藏 布 
江 一 线 大 部 、 昌 都 市 西南 部 以 及 隆子 站 的 干旱 致 灾 
危险 性 最 高 ;林芝 市 中 部 和 聂 拉 木 站 、 错 那 站 处 于 
中 等 危险 性 ,其 他 地 区 属于 较 高 等 级 。 抽 种 -成 熟 
期 , 除 雅 鲁 藏 布 江 一 线 大 部 和 错 那 站 、 隆 子 站 、 八 窒 
站 外 ,其 他 地 区 旱灾 的 危险 性 多 属于 中 等 级 别 以 


4y Be ph oh — Beta Ht 37.9 37.9 
of fay at 13.9 51.8 
第 三 载荷 量 8.1 59.9 
第 四 载荷 量 6.6 66.5 
第 五 载荷 量 5.6 72.1 
全 生育 期 第 一 载荷 量 35.8 35.8 
SP at 24.4 60.2 
第 三 载荷 量 7.3 67.4 
第 四 载 和 荷 量 7.1 74.6 
第 五 载荷 量 4.7 79.3 


下 ,其 中 胎 拉 木 站 和 芒 康 站 危险 性 最 低 。 全 生育 
期 ,干旱 致 灾 人 危险 性 等 级 均 处 于 中 等 以 下 ,林芝 站 
和 米 林 站 最 低 ; 中 等 危险 性 区 域 主要 位 于 雅 鲁 藏 布 
江 一 线 , 其 他 区 域 危险 性 较 低 。 


3 讨论 


(1) 结合 前 面 的 干旱 过 程 覆 盖 范 围 ,分 析 各 种 
空间 模 态 干旱 强度 随时 间 的 变化 特征 。 播 种 -分 划 
期 ,EOF1 在 1990.1992 .1997 年 和 2015 年 这 4a 均 为 
干旱 强 年 。 分 蓝 - 抽 穗 期 , EOF1 在 1982 1986— 
1987 .2000 年 和 2008 一 2009 年 其 时 间 系 数 绝对 值 较 
高 ,结合 图 4b ,发 现 1982 1986 一 1987 年 和 2009 年 
为 明显 的 干旱 强 年 ,其 中 2009 年 的 时 间 系 数 绝 对 值 
最 高 ,当年 对 应 的 干旱 覆盖 范围 也 大 (84% ) ,这 与 袁 
雷 等 下 对 西藏 西部 地 区 雨季 得 出 的 “干旱 强度 在 
1981 一 1990 年 和 2001 一 2010 年 偏 高 ,在 2009 年 达 
到 峰值 ”的 结 轮 较 一 致 ;EOF2 在 1993、1995 年 和 
2018 年 为 干旱 强 和 年。 抽穗- 成熟 期 ,结合 图 4c , 发现 
EOF1 7E 1983 , 1986 , 1992 4£ FI 2009 fF 3X 4 a 为 突出 
的 干旱 强 年 ,其 中 1983 年 的 时 间 系 数 绝 对 值 最 高 ， 
同时 干旱 影响 范围 为 最 大 值 (84%);EOF2 在 1982 一 
1983、1986 年 和 1992 年 为 典型 的 干旱 强 年 ; 近 40a 
来 ,抽穗 -成 熟 期 2 个 模 态 的 干旱 强 年 与 袁 雷 等 ”对 
西藏 沿江 一 线 雨 季 得 出 的 “干旱 强度 在 1981—1990 
年 达到 最 大 ”的 结 轮 较 一 致 。 全 生育 期 ,结合 图 4d， 
发 现 EOF1 在 1982—1983 , 1986—1987 , 1992 年 和 
2009 年 为 显著 的 干旱 强 年 。 

(2) 文中 增加 了 2014 一 2020 年 时 间 序 列 ,得 出 
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图 例 EOF1 

。 气象 站 点 号 <-0.1 

DRF = -0.1~0.0 m >0.2 
= 0.0~0.1 


ES || 


(c) 4d 58 THRE SHEOF2 


图 例 EOF2 
。 气 象 站 点 <01 mm 0.1~0.2 0 180 km 
ORR = —0.1~0.0 m >0.2 
= 0.0~0.1 


(e) 抽穗 -成 熟 期 EOF2 yN 


图 例 EOF2 


ep <0.l m0.1~0.2 0 180 km 
"和 气象 站 后 = —0.1~0.0 m >0.2 — 
CI 


= 0.0~0.1 
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(b) 分 药 - 抽 穗 期 BOF1 人 


图 例 EOF1 
。 气象 站 点 <0.1 
口上 县 界 


== 一 0.1~0.0 >0.2 L 


= 0.0~0.1 


(d) 抽穗 -成 熟 期 EOF1 


图 例 EOF1 ð 

. fyb 回 <-0.1 m0.1~0.2 

a =-0.1~0.0 ma>0.2 
mm 0.0~0.1 


(全 生育 期 EOF1 


EOF1 
图 例 <-0.1 mm 0.1~0.2 0 


k= fp,» C4 
。 气 象 站 点 = —0.1~0.0 ma >0.2 
口上 县 界 ma 0.0~0.1 


180 km 


注 :EOF1、EOF2 分 别 为 经 验 正 交 函 数 法 (EOF ) 的 第 一 模 态 .第 二 模 态 。 下 同 。 


图 7 1981 一 2020 年 西藏 青 称 主 要 种 植 区 各 生育 期 干旱 强度 EOF 不 同 模 态 空间 分 布 特征 
Fig.7 Spatial distribution characteristics of different modes of drought intensity EOF in different growth stages 


of highland barley growing areas in Tibet from 1981 to 2020 


1981—1990 F TFA tet A] ee .持续 日 数 最 长 、 
发 生 频 率 最 大 ,2001 一 2010 年 均 最 小 ( 短 ) , BRE LR 
楠 " 等 得 出 “旱灾 在 1991 一 2005 年 发 生 较 为 频繁 ， 
2002 年 后 每 年 发 生 频 次 均 低 于 多 年 平均 值 ” 的 结论 
不 尽 相 同 ,分 析 发 现时 间 序 列 的 增加 并 未 改变 干旱 
频率 的 低 发 时 间 段 ,出 现 高 发 频率 干旱 年 份 不 一 致 
的 原因 可 能 是 本 文 研究 仅 涉 及 全 年 的 部 分 月 份 
(4 一 9 月 ) ,虽然 1 一 3 月 和 10 一 12 月 的 和 干旱 灾害 发 
生 相 对 较 少 ,但 也 会 发 生 。 不 同 生育 期 干旱 发 生 频 
率 的 高 值 区 主要 集中 在 日 喀 则 市 东部 .拉萨 市 南 
部 、 昌 都 市 西南 部 和 山南 市 中 部 地 区 , 低 值 区 位 于 


林芝 市 中 部 ,这 与 能 俊 楠 等 "得 出 的 区 域 范围 较为 
吻合 ,但 与 袁 雷 等 ”对 西藏 雨季 干旱 频率 的 结论 有 
所 差别 ,原因 可 能 是 时 间 序 列 \ 人 研究 区 域 .干旱 指标 
的 选取 不 同 造成 。 

(3) 综合 4 个 生育 期 干旱 致 灾 危 险 性 等 级 , 西 
藏青 牧 主 要 种 植 区 干旱 致 灾 危 险 性 较 蜗 的 区 域 主 
要 分 布 在 日 喀 则 市 中 东部 、 拉 醒 市 南部 和 山南 市 北 
部 ,该 区 域 多 位 于 雅鲁藏布江 一 线 ,降水 相对 偏 少 ， 
蒸 散 较 大 ,更易 发 生 干 量 ; 较 低 区 域 主要 集中 在 林 
艺 市 ,该 区 域 属 于 热带 湿润 和 半 湿 润 气候 ,降水 丰 
富 、 素 发 量 相对 较 小 ,干旱 事件 相对 少 发 。 熊 俊 楠 
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时 间 系 数 


时 间 系 数 


(a) 播种 -分 药 期 Rn eh (b) 488 Hi 


图 例 致 灾 危 险 性 等 级 一 图 例 致 灾 危 险 性 等 级 ” 
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图 8 1981 一 2020 年 西藏 青 牺 主 要 种 植 区 各 生育 期 干旱 强度 EOF 不同 模 态 对 应 的 时 间 系 数 
Fig. 8 Time coefficient of EOF of different modes in different growth stages of highland 
barley growing areas in Tibet from 1981 to 2020 


ez 


。 气 象 站 点 回 <10 mm 20-25 Be 。 气象 站 点 <10 mm 20~25 0 180km 
ORR m1015 m>25 ooo o 县 界 = 10~15 mm >25 
== 15~20 = 15~20 


(c) 抽穗 -成 熟 期 (d 全 生育 期 


Zz 


图 例 致 灾 危 险 性 等 级 


El fil 致 灾 和 危险 性 等 级 
。 所 加 <10 mm20~25 0 180km 。 气 象 站 点 <30 60-75 0 L 

He S 国 10-15 >25 一 器 AR 国 30-45 >75 
ORF 15-20 mm 45~60 


图 9 1981—2020 F PRA R ERPE KK EA WT BIC I Bn EG BIZ TH] PA 
Fig. 9 Spatial distributions of drought hazard grades in different growth stages of highland 
barley growing areas in Tibet from 1981 to 2020 
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等 "基于 历史 气象 灾害 事件 、 史 继 清 等 "利用 自然 
水 分 亏 缺 率 指数 对 西藏 地 区 气象 干旱 空间 分 布 特 
征 所 得 出 的 结论 与 本 文 结 论 较 一 致 ; 而 金 建 新 等 ” 
的 研究 结 采 认为 “干旱 指数 从 西 问 东 逐 渐 递 增 , 从 
西部 地 区 的 特 干 旱 气 候 区 依次 向 东部 的 湿润 半 湿 
润 及 湿润 气候 区 过 渡 ” ,与 本 文 结论 有 所 差异 , 究 其 
原因 在 于 选取 的 指标 不 一 致 .两 者 的 计算 方法 存在 
差异 等 。 研 究 区 分 井 - 抽 穗 期 致 灾 危 险 性 总 体 较 
高 ,可 能 由 于 分 碍 -抽穗 期 是 至 穆 需 水 关键 期 , 且 此 
WAKER REER 气温 回升 快 ,最终 导 致 
了 该 生育 期 干旱 最 为 严重 ,进而 证 实 了 熊 俊 楠 等 "得 
出 “旱灾 多 发 生 在 上 半年 ,从 1 月 开始 逐渐 递增 ,在 6 
月 达到 峰值 "的 研究 结论 。 


4 结论 


(1) 1981 一 2020 年 西藏 青 穆 主 要 种 植 区 的 气温 
呈 显 著 上 升 趋势 ,降水 呈 不 显著 增加 趋势 ,季节 性 
气温 和 降水 变化 趋势 差异 明显 。 多 年 来 秋冬 季 降 
水 呈 减 少 趋势 , 春 夏 季 则 呈现 不 同 速率 的 增加 趋 
势 ; 四 季 的 年 际 气 温 变化 均 呈 显著 上 升 趋势 ,其 中 
冬季 气温 上 升 速率 最 快 ,达到 0.46 C+ (10a). 

(2) 1981 一 2020 年 西藏 青 称 主要 种 植 区 不 同 生 
育 期 干旱 覆盖 范围 和 持续 日 数 总 体 呈 现 不 同 程度 
的 下 降 趋势 ,干旱 强度 均 有 所 增强 ( 除 播 种 -分 幕 期 
外 ) ;分 葛 - 抽 穗 期 干旱 频率 最 高 ,播种 -分 幕 期 次 
之 ,抽穗 -成 熟 期 最 低 。 

(3) 通过 EOF 分 析 的 西藏 青 穆 主 要 种 植 区 干旱 
强度 分 布 特点 :播种 -分 雍 期 呈现 中 部 与 中 西部 相 
反 ;分 昔 - 抽 穗 期 呈现 全 区 一 臻 .中 部 与 中 西部 相 
反 ; 抽 穗 -成 熟 期 呈现 全 区 一 致 ( 除 野 拉 木 站 外 ) \ 西 
机 与 东北 相反 ;全 生育 期 呈现 全 区 一 致 的 多 种 模 
态 。 中 部 与 中 西部 相反 的 典型 年 份 为 2015 年 ( 播 
种 -分 碍 期 , 定 日 站 、 拉 和 孜 站 严重 干旱 ) 和 1995 年 (分 
幕 -抽穗 期 ,昌都 市 大 部 、 林 艺 市 东部 严重 干旱 ) , 西 
南 与 东北 相反 典型 年 份 为 1986 年 (日 喀 则 市 东部 严 
Bt). 

(4) Pe HER FS PRE EPEK FR BK SE EG 
ERR-I BET Aad BEF EY, ES Ee SL ah 
低 两 边 高 的 格局 特征 ;处 于 抽穗 -成 熟 期 和 全 生育 
期 ,整体 上 呈现 西高 东 低 的 的 分 布 特征 。 
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Spatiotemporal distribution of drought and hazard assessment 
of highland barley in Tibet 


SHI Jiqing', GAN Chenlong, ZHOU Kanshe', YUAN Lei’, ZHANG Dongdong’ 
(1. Tibet Climate Center, Lhasa 850000, Tibet, China; 2. Maizhokunggar Meteological Bureau, Lhasa 850000, Tibet, China) 


Abstract: In Tibet, highland barley is mostly cultivated on the primary crop-producing area in Tibet, China. How- 
ever, drought resulting from climate change and intense human activities agriculture is considerably affecting ag- 
riculture in this area. The daily meteorological drought comprehensive index (MCI) of the growth period of 25 
meteorological stations in the main growing areas of highland barley in Tibet, was calculated according to the na- 
tional standard of “Meteorological Drought Grade” (GB/T 20481-2017), and the spatial and temporal variation 
characteristics of meteorological drought and the distribution of disaster risk were analyzed in the main growing 
areas of highland barley in Tibet. The results revealed that: (1) In the last 40 years, the air temperature in the main 
growing areas of highland barley in Tibet increased significantly (P<0.01), whereas the precipitation did not in- 
crease significantly. The seasonal variation trend of temperature and precipitation differed considerably. (2) The 
coverage and duration of drought exhibited a decreasing trend of various degrees, and the drought intensity (ex- 
cept sowing-tillering stage) enhanced. Drought was the most severe at the tillering-heading stage and in 1981— 
1990, whereas drought was the least severe in 2001—2010. (3) The high-value areas of drought frequency in vari- 
ous growth stages were in eastern Xigaze, southern Lhasa, southwestern Qamdo, and central Shannan, whereas 
low-value areas were in central Nyingchi. (4) The patterns of empirical orthogonal function analysis in various 
growth stages were opposite in central and central-western regions and consistent in the whole region, and oppo- 
site in southwest and northeast regions. (5) The drought disaster risk grades of each growth period were higher in 
the west and lower in the east. The disaster risk grades of the tillering-heading stage were generally higher, that of 
the heading-maturity stage was lower, and that of the sowing-tillering stage was the lowest. 


Key words: MCI index; hazard risk; drought intensity; barley; Tibet 


